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INTERPOLATION FOR USE IN CHANNEL ESTIMATION 



Technical Field 

The present invention is related to telecommunications 
5 systems and in particular it concerns a method and a device 
for the estimation of the transfer function of a transmission 
channel . 
Background Art 



10 mobile radio communication systems is the W-CDMA (Wideband - 
Code Division Multiple Access) technique, by virtue of its 
high spectral efficiency as compared to other multiple access 
techniques. Within said technique and in particular in the 
FDD (Frequency Division Duplexing) mode different solutions 

15 have been adopted to increase system capacity* Among said 
solutions it is worth mentioning the coherent tracking of the 
signal received by the radio mobile station or mobile 
terminal . 



20 estimation of the transfer function of the transmission 
channel, referenced to in the sequel as "channel estimation" , 
in the presence of fading and multiple reflections due to 
user's motion in the covering area, without giving any 
details about a user's displacement velocity. 

25 Like other system functions, also the channel estimation 

function may be assigned to processing devices (for example, 
DSPs, microcontrollers, etc) and in such a case this is 
achieved by means of software or through dedicated devices 
(for example FPGAs, ASICs, etc.) and therefore the 

30 implementation modality is mainly hardware based. During the 
definition of the telecommunications system architecture, the 
design engineer has to find out the best distribution of the 
different tasks among programmed devices and wired devices. 



As is well-known, one of the most used techniques in the 



This type of signal tracking requires a very accurate 
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Thus, the different architectures achievable must then 
be analysed in terms of computational burden for a DSP 
device, traffic volume over the bus dedicated to data 
transfer, and the optimum distribution among hardware and 
5 software resources . 

Considering by way of an example a complex system such 
as a modem for a UMTS base station or mobile, terminal., it is - 
evident that the communication bus between the DSP device and 
the hardware resources acts as a bottleneck of the system. 
10 For instance, as to a radio base station, the load of the bus 
becomes greater and greater as the number of users increases, 
and when the functions involved require the transfer between 
DSP and hardware resources of considerable bulks of data. 

Channel estimation is one of said functionalities, since 
15 it has to be repeated for any fingers (sub- channels) of each 
individual user . 

Let us consider the worst case in which the maximum 
number of users Nu = 128 are communicating within a cell and 
each user has a finger number Nf = 8 for each Rake receiver 
20 (a receiver which is typical for example of UMTS base 
stations) . The quantity of data to be transferred for each 
finger in a time slot (time interval) of duration Ts = 666 ms 
may be estimated to be equal to a number of bits NB = 400. As 
a consequence, there will be a considerable data stream over 
25 the bus, namely: 

Nb-Nf *Nu/Ts=400-8-128/666aO-6 = 615 Mbit/s 
As of today, various algorithms capable of performing 
the channel estimation are known; among them there are for 
instance those described in the document EP 0 912 019, 
30 wherein two different interpolation methods, the one of 
linear type, the other through Kalman filters, are used, the 
method described in the document EP 1 032 168, wherein the 
interpolation is performed with the method of Lagrange's 
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polynomials, and that described in the document US. .5 , 8 86 , 911 , 
in which use is made of the bi- sect ion method. 

The solutions described in the above documents have 
however some drawbacks and are not directly applicable to 
5 UMTS systems, such as the bi- section method, or are extremely 
complex for an implementation of hardware type. 

Taking into account as a matter of fact, the 
considerable data stream already present on the bus of a base 
station, it should be preferred not to further overload such 

10 a bus, since this would require complex interpolation 
operations within the estimation of the channel, such as the 
operations necessary to perform an interpolation by the 
method of Lagrange's polynomials. 

The device can be fully implemented through hardware and 

15 therefore can be easily integrated into a Rake receiver for 
base stations, keeping the same performance as of the 
software solutions based on DSP processors. Since no data 
transfer with the DSP processor is required, the 
communication bus is left free for other tasks. 

20 A particular subject matter of the present invention are 

a method and a device for the estimation of the transfer 
function of a transmission channel, as described in the 
appended claims . 
Summary of the invention 

25 The method and the device subject matter of this 

invention make use of an algorithm of a low complexity, for 
the estimation of the transfer function of a transmission 
channel, suitable for both the transmission path toward a 
base station, called "up-link", and for the transmission path 

30 toward a mobile terminal, called " down- 1 ink " . 

The method according to the invention envisages the 
computation of a plurality of channel coefficients, included 
among known channel coefficients corresponding to pilot 
symbols, through the reiteration of an interpolation 
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algorithm / capable of calculating an intermediate point 
between a first extreme and a second extreme of a determined 
interval. The first extreme is formed by at least two known 
points and the second extreme is formed by at least one known 
5 point • 

The intermediate point to be calculated has an abscissa 
equal to the abscissa value of_ a mean point between the 
points defining the interval rounded off to the integer 
closest to the first extreme, and has an ordinate equal to 
10 the arithmetic average between the ordinate of the known 
point of the second extreme and the ordinate of a point, 
chosen between the two known points of the first extreme, 
having a distance from the intermediate point equal to the 
distance between the intermediate point and the known point 
15 of the second extreme. 

Brief Description of Drawings 

This and other characteristics of this invention will 
become evident from the following description of a preferred 
embodiment of the same, given by way of a non- limiting 
20 exanple, and from the attached drawings, wherein: 

Figure 1 is a schematic representation of the symbols of 
a data channel DPDCH and of a control channel DPCCH, in the 
case of a transmission path toward a radio mobile station; 

Figure 2 is a graph illustrating a plurality of channel 
25 coefficients obtained through linear interpolation between 
two known channel coefficients. 

Figure 3 is a flow chart illustrating a first 
interpolation algorithm implemented according to this 
invention; 

30 Figure 4 is a schematic representation illustrating a 

sequence of operations for computing the channel coefficients 
through the algorithm of Figure 3, when N PJLOT = 3; 
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Figure 5 is a schematic representation illustrating a 
sequence of operations for computing the channel coefficients 
through the algorithm of Figure 3, when Npilot - 4; 

Figure 6 is a schematic representation illustrating a 
5 sequence of operations for computing the channel coefficients 
through the algorithm of Figure 3, when N PIK>T = 5; 

Figure 7 is a schematic representation illustrating a 
sequence of operations for computing the channel coefficients 
through the algorithm of Figure 3, when Nph^t = 6; 
10 Figure 8 is a schematic representation illustrating a 

sequence of operations for computing the channel coefficients 
through the algorithm of Figure 3, when N PILOT = 7; 

Figure 9 is a schematic representation illustrating a 
sequence of operations for computing the channel coefficients 
15 through the algorithm of Figure 3, when N PILOT = 8; 

Figure 10 is a state diagram summarising the operations 
required for computing the channel coefficients in the cases 
depicted in Figures 4 to 9; 

Figure 11 is a schematic representation illustrating a 
20 sequence of operations for computing the channel coefficients 
through a second algorithm implemented according to the 
present invention, when N PIIlCyr = 6; and 

Figure 12 is a block diagram of a circuit performing an 
interpolation algorithm, implemented according to the present 
25 invention. 

Detailed description of the preferred embodiments 

The method and the device according to the present 
invention will now be described in detail with reference to 
the UMTS (Universal Mobile Telecommunications System) 
30 systems, relating to the up-link path in the FDD mode. 

The coherent tracking for the radio interface of UMTS 
systems may be obtained by using time -multiplexed pilot 
symbols in the control channel DPCCH (Dedicated Physical 
Control Channel) . 
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In Figure 1 there is schematically depicted the time 
relation between the symbols of the data channel DPDCH 
(Dedicated Physical Data Control Channel) and the symbols of 
control channel DPCCH in a time interval called "time slot" . 
5 In the case illustrated in said figure, there are at the 
beginning of the slot three pilot symbols, Xq(0) , Xq(1) and Xq 
(2) , followed by other symbols used _ for the. control, in 
particular four symbols TFCI (Transport - Format Combination 
Indicator) , two FBI (Feedback Information) symbols and a TPC 

10 (Transmit Power Control) .symbol. 

Pilot symbols XqOO , with k«0,l,2,.., N PILO t-i (where Npilot 
is the number of pilot symbols in the slot) are also known to 
the receiver which utilises such information jointly with the 
received symbols Z(k) f in order to evaluate the channel 

15 effect on each pilot symbol Xg(k) . 

By carrying out the complex product of the received 
symbol Z(k) by the corresponding pilot symbol Xg(k) one 
obtains the channel coefficients, in particular the in phase 
and quadrature components of the product correspond to the 

20 components C Q e C x of the channel coefficients 
C(k)=c I (k)+jc Q (k) . 

It becomes therefore evident that it is possible to 
estimate the channel coefficients only in relation of the 
pilot symbols transmitted over channel DPCCH. However, since 

25 the pilot symbols are not transmitted continuously but only 
in coincidence with the initial part of a time slot, 
performing the channel estimation requires the computation, 
through some form of interpolation, of the channel 
coefficients corresponding to the missing pilot symbols. 

30 The number of pilot symbols within an individual slot 

may vary from a minimum of three to a maximum of eight; 
therefore the number of channel coefficients that will have 
to be estimated through interpolation, shall vary between 
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seven and two, being ten the total number of symbols present 
in a slot of a control channel DPCCH. 

Figure 2 depicts for instance the case in which seven 
channel coefficients (in the figure are represented the in 
5 the phase components only) , globally indicated with reference 
2, have been computed through interpolation between the last 
coefficient Ci (Npn/jr-l) of the current slot, corresponding to 
a value of abscissa value A and a first coefficient of the 
following slot, indicated in the figure as C x (10) and 

10 corresponding to a value of abscissa B. The pilot symbols , 
thus also the channel coefficients, always correspond by 
common assumption to integer abscissa values {0, 1, 2, 3, 
,..) on a horizontal time axis (Time)* 

The interpolation method and device implemented 

15 according to the present invention allow computation of the 
components C;c<k) e C Q (k) of the channel coefficients without 
using complex operations such as multiplication and division, 
but only utilising additions and divisions by two (the latter 
ones of easy implementation in hardware through a shift 

20 operation on the right side of a register) allowing a 
considerable reduction in the hardware architecture 
complexity of the interpolation unit. 

The interpolation method according to the present 
invention will now be described in detail with reference to 

25 the flow chart of Figure 3 , The algorithm shown in Figure 3 
allows the computation, through interpolation of a plurality 
of points, included between a last channel coefficient of 
abscissa A, a current slot L, and a first channel coefficient 
of abscissa B, of a slot L+l subsequent to said current slot. 

30 Since the minimum number of channel coefficients in each slot 
is equal to three, the coefficient, of abscissa A-l, 
immediately preceding the last channel coefficient of each 
slot is always known and therefore can be used for the 
interpolation computation. 
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Within the algorithm of Figure 3 some variables^ X, Y, 
Z, W are used with the following meaning: 

X: abscissa of the known left-hand point of the 
interpolation interval, a value initially corresponding to A; 
5 Y : abscissa of the known right-hand point, used for the 

interpolation, a value initially corresponding to B; 

Z: abscissa of an intermediate point between X and Y; 

W: abscissa of the extreme left point, actually used for 
the interpolation, corresponding, depending on the case, to 
10 point X or to point X-l; 

f (kj : ordinate corresponding to a generic abscissa k. 

Variables X e Y are initialised as being equal to the 
extremes A and B in the initial block indicated by reference 
32 in Figure 3 . 

15 Subsequently in block 36 of such figure, abscissa Z is 

obviously calculated as an integer of the intermediate point 
between X an Y. The value of abscissa Z is rounded off to the 
lower integer through function FLOOR (symbolically indicated 
in the figure) . In the same block the computation is made of 

20 the abscissa of the extreme left-hand point W, which will be 
equal either to X or to X-l, depending on whether the value 
of X+Y is odd or even, respectively (by applying the formula 
W=2«Z-Y) . Then the ordinate f(Z) of the intermediate point Z 
is calculated by arithmetically averaging between the 

25 ordinate of the extreme left-hand point W and the ordinate of 
the right-hand extreme Y, which are known. 

Then, still within block 36, a check is made on whether 
the just calculated point on abscissa Z corresponds to the 
abscissa point X + 1, If such a condition is not met, value Z 

30 is assigned to abscissa Y of the known right-hand point, and 
again operations contained in block 36 are applied. In 
practice Block 36 algorithm is recursively applied to the 
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half -interval at the left side of the point previously 
computed, until the abscissa point X+l is attained. 

Subsequently, one should consider the operations 
indicated as a whole with reference 38, operating by 
5 increasing abscissas, starting from the last computed point 
X+l, and searching a first point still to be calculated. Once 
such a point has been found, the operations contained in 
block 4 0 are carried out or one terminates the algorithm if 
the right-hand extreme B of the interpolation interval has 
10 already been reached, for all the points have already been 
computed. 

If however a point has been found which is still to be 
computed and whose abscissa is equal to the value taken by 
variable 2, the next step is to assign to variable X the 

15 value Z-l that represents the first known left-hand point, 
and to search, by subsequent attempts, for the first known 
right-hand point of Z, that is. the value which is assigned 
to variable Y. To perform the latter operation, variable Y is 
increased by a unit each time, until a corresponding known 

20 coefficient has been reached. These operations are carried 
out within block 40 

With these new values of X and Y, the recursive 
algorithm of block 36 is again applied, and in cascade use is 
made of the procedure contained within block 38 for the 

25 search of a new point to be computed. If all the points have 
been computed, the algorithm terminates; otherwise a new X-Y 
interval is determined to which again recursive algorithm 36 
is to apply. 

The previously-mentioned algorithm, shown in Figure 3, 
30 allows computation, exclusively by sums, divisions by two and 
compare operations between registers, of any number of 
channel coefficients contained in an interval defined on the 
left side by two known coefficients (of abscissa A and A-l) 
and on the right side by a known coefficient (of abscissa B) . 
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For a better understanding of the general interpolation 
algorithm described above, we will now consider its 
application to six particular cases, typical of the UMTS 
system, from the case in which the pilot symbols (Npilot) are 
5 three, and the coefficients to be computed are seven, up to 
the case where the pilot symbols are eight and the channel 
coefficients to be computed are two only. _ _. . 

The fist case (Np ILOT =3) is schematically illustrated in 
figure 4. On the abscissa the channel coefficients are 

10 plotted corresponding to the three pilot symbols of the 
current slot, namely to abscissa positions 0, 1 (A-l) and 2 
(A), and. to the first pilot symbol of the subsequent slot, to 
the abscissa position 10 (B) respectively. The channel 
coefficients to be computed are therefore those corresponding 

15 to the abscissa positions 3 to 9 . 

The interpolation operation is carried out by subsequent 
steps; the number of steps directly depends upon the number 
of channel coefficients to be computed. In case of Figure 4 
the computation is performed in seven steps (Step 1 to Step 

20 7). 

At Step 1, the algorithm of block 36 of flow chart of 
Figure 3 is applied by assigning to variables X and Y the 
values A (which is worth 2) and B (which is worth 10) , 
respectively: 
25 X=A=2 

Y=B=10 

Thus the following is obtained: 

Z=FLOOR [ (2+10) /2] =6 ; abscissa of the intermediate 
point to be computed (denoted by letter C in the Figure) . 
30 W=2-6-10=2 ; abscissa of the left-hand extreme' point, 

corresponding in this case to X. 

Once the abscissa of the extreme left-hand point W, of 
the extreme right-hand point Y and of the intermediate point 
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Z, of which the ordinate is to be computed, are known, the 
next operation is to calculate the value of the corresponding 
coefficient C (6) =Ci (6) + jC Q (6) f computing the components C x e 
C Q by the arithmetic averaging rule (point C on the graph of 
5 Figure 4) : 

Cj(6) = [C I (2)+Ci(10) ] /2 and C Q (6) = [C Q (2) +C Q (10) ] /2 
After checking that abscissa Z of the computed point , in 
this case being equal to 6, still does not correspond to the 
point X+l=3, one assigns the value Z=6 to variable Y and 
10 determine a new coefficient (point D, step 2) : 
X=2; Y=6 

Z=FLOOR [(2+6)/2]=4 
W=2-4-6=2 

Ci (4) = [Ci (2) +d (6) ] /2 and C Q (4) = [C Q (2 ) +C Q (6) ] /2 
15 Also in this case the output condition from block 3 6 is 

not met, since Z=4 is different from X + 1= 3; thus the next 
step is to compute a new coefficient (point E, step 3) : 

X=2; Y=Z=4 

Z=FLOOR [ (2+4) /2] =3 
20 W=2-3-4=2 

C z (3) - [C t (2) +d (4) ] /2 and C Q (3) = [Cq (2 ) +C Q (4) ] /2 

At this point a verification is made of condition Z=X+1 
which allows passing from block 36 to block 38. 

By applying the rules described within block 3 8 , one 
25 finds the first point still to be computed on the right hand 
of the point of abscissa Z=3 is the abscissa point 5, then 
one realises that said point does correspond not to the right 
hand extreme B=10 and go on to block 40, having assigned 
value 5 to the variable Z. 
30 Within block 4 0 one determines the first known left-hand 

point X=Z-1=4 and the first known right-hand Y=6, and with 
these values of X and Y one then goes back to block r 36 to 
calculate a new coefficient (point F, step 4) : 
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X=4; Y=6 

Z=FLOOR [(4+6)/2)=5 
W«2-5-6=4 

Ci (5) = [C x (4) +Ci (6) ] /2 and C Q (5) = [C Q (4 ) +C Q (6) ] /2 
5 In this case the condition Z=X+1 is immediately 

verified, and one goes on again to block 3B where the first 
right-hand point still to be computed -is the abscissa point 
Z=7. Within block 40, the first known left-hand point X=Z-1=6 
and the first known right-hand point Y=10 are determined, and 
10 with these values of X e Y one goes back to block 36, to 
compute a new coefficient (point G, step 5) : 
X-6; Y=10 

Z=FLOOR [ (6+10) /2] =8 
W=2-8-10=6 

Ci(8) = [Ci(6)+Ci(10) ]/2 and C Q (8 ) = [C Q (6) +C Q (10) ] /2 
Condition Z= X+l is not verified; thus one remains 

within block 3 6 to calculate a new coefficient (point H, step 

6) : 

X=6; Y=Z=8 
20 Z-FLOOR [(6+8)/2]=7 

W=2-7-8 = 6 

Ci (7) - Id (6) +Cx (8) ] /2 and C Q (7) « [C Q (6) +Cq(8) ] /2 
Having reached point Z=X+1*7, one goes on to block 38 
where it must be verified whether the right -end extreme has 
25 not yet been reached, and a determination of abscissa Z = 9 
of the next point to be calculated is to be effected. Within 
block 40 the first known left-hand point X=Z-1=8 and the 
first known right-hand point Y=10 will then be determined, 
and with these values of X and Y we go back to block 36 to 
30 compute a new coefficient (point I, step 7) . 
X=8; Y=10 

Z=FLOOR [ (8+10) /2] =9 
W~2-9-10«8 
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Ci (9) = [C I (8)+Ci (10)]/2 and G Q (9) - [C Q (8) +C Q (10) ] /2 
Condition Z=X+1 has now been verified; therefore we go 
on to block 3 8 where variable X is increased to reach the 
value Z=10 which, as corresponding to the right-hand extreme 
5 B, leads to termination of the interpolation algorithm (step 
"STOP" in Figure 3) . 

As may be derived from the above equations, all the, 
channel coefficients have been calculated only by sums and 
divisions by two, i.e, operations which are easily carried 
10 out directly on hardware* 

In the case (Npiiot-3 ) . described above, the known 
coefficient of abscissa A-l has not been used, since the 
distance between points A and B is a power of 2, in 
particular 2 3 =8. This condition is brought about also in the 
15 case Npi L 6 T =7 / where such distance is 2 2 =4 . 

The algorithm depicted in figure 3 allows however the 
computation of all the intermediate coefficients also in 
cases when such a distance is not a power of 2 , as in the 
remaining cases, illustrated in the sequel. 
20 In Figure 5 there is depicted the case in which the 

pilot symbols are four (Npi LO t=4) and the coefficients to be 
calculated are six, from C(4) to C(9). The interpolation plan 
illustrated in figure 5, as will be appear from the following 
detailed explanation, uses the abscissa point A-l as the 
25 left-hand extreme for the interpolation. Even if it 
introduces a slight degradation of the overall 
characteristics of the system, this approximation greatly 
simplifies the computation of the interpolated coefficients. 

Let us now see how the algorithm of figure 3 is applied 
30 to the present case (Npilot^) . 

As the first step, use is made the first time of 
algorithm of block 3 6 of Figure 3 flow chart, by assigning 
the values of A (which is worth 3) and of B (which worth 10), 
to variables X and Y, respectively: 
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X=A (A=3) 
Y=B . (B=10) 
Thus : 

Z=FLOOR [ (3 + 10) /2] «6; abscissa of the intermediate point 
5 to be calculated (denoted in the figure by letter C) ; 

W=2-6-10=2; abscissa of the extreme left-hand point, in 
this case it does not correspond, to X. 

Once the abscissas are known of the extreme left-hand 
point W, the extreme right-hand point Y and the intermediate 
10 point Z of which the ordinate is to be calculated, one starts 
computing the value of the corresponding coefficient 
C(6) =C!(6)+jC Q (6) , calculating components C z e C Q by the 
arithmetic averaging rule (point C in the graph of Figure 5) : 
C x (6) - [Cj (2) +d (10) ] /2 and C Q (6) = [C Q (2) +C Q (10) ] /2 
15 After verifying that abscissa Z of the calculated point , 

in this case equal to 6, does not still correspond to point 
X+l=4, one assigns the value Z=6 to variable Y, and a new 
coefficient (point D, step 2) is computed: 
X=3; Y«6 
20 Z=FLOOR [(3+6)/2]*4 

W=2-4-6=2 

C I (4)«[C I (2)+Ci(6)]/2 and Cq<4) = [Cq(2) +Cq(6) ] /2 
In this case the output condition of block 36 is 
verified, since Z=4 is equal to X+l=4 ; consequently one goes 
25 on to block 38. 

By applying the rules described within block 38, one 
finds that the first point still to be calculated on the 
right hand of abscissa Z=4, is the abscissa point 5, verifies 
that such point does not correspond to the right-hand extreme 
30 B=10, and goes on within block 40 having assigned the value 5 
to variable Z, 

Within block 40, one will determine the first known 
left-hand X=Z-1=4 and the first known right hand point Y=6, 
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and with these values of X and Y around" block 3 6 where a new 
coefficient (point E, step 3) is computed. 
X=4; Y=6 

Z= FLOOR [(4 + 6)/2]=5 
5 w=2-5-6=4 

Ci (5) = [C x (4) +d (6) ] /2 and C Q (5) = [C Q (4) +C Q (6) ] /2 
At this point the condition Z=X+1 must be verified that 
allows passing from block 36 to block 38. 

By applying the rules described within block 38, one 
10 finds that the first point still to be calculated on the 
right-hand of abscissa point' Z=5 is the abscissa point 7, 
verifies that such point does not correspond to the right- 
hand extreme B=10, and goes on within block 40, having 
assigned the value 7 to variable Z. 
15 Within block 40 , the first known left-hand point X=Z-1=6 

and the first known right-hand point Y=10 are determined, and 
with these values of X and Y one returns to block 36 to 
calculate a new coefficient (point F, step 4) : 
X=6; Y=10 
20 Z=FLOOR [ (6+10) /2] =8 

W=2-8-10 = 6 

Ci(8) = fCi(6) +Ci(10) ] /2 and C Q (8) - [C Q (6) +C Q (10) ] /2 
Condition Z=X+1 is not verified, therefore one remains 
within block 3 6 to calculate a new coefficient (point G, step 
25 5) : 

X=6; Y=Z=8 

Z=FLOOR [(6+8)/2]=7 

W=2-7-8=6 

Ci (7) = [C x (6) +C X (8)3/2 and C Q (7) = [C Q (6) +C Q (8) ] /2 
30 In this case condition Z=X+1 is immediately verified, 

and one returns to block 3 8 and finds out that the first 
right-hand point still to be calculated is the abscissa point 
Z=9. Within block 40 the first known left-hand point X=Z-1=8 
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and the first known right-hand point Y= 10 are determined, and 
with these values of X and Y one returns to block 36 to 
compute a new coefficient (point H, step 6) ; 

X=8; Y=10 
5 Z=FLO0R [(8+10)/2]=9 

W=2-9-10 = 8 

_ _ Ci (9) - [Cr (8) +d (10) J/2 - and C Q (9) = [Cg(8) +Cq (-10) ]-/2 

Condition Z=X+1 is now verified; therefore one goes on 
to block 3 8 where variable Z is increased passing to value 
10 Z=10, which, corresponding to the right-hand extreme Y, leads 
to terminate the interpolation algorithm (step STOP in Figure 
3) . 

Figure 6 illustrates the case where each slot contains 
five pilot symbols (Npilqt-5) , and the number of channel 

IS coefficients to be calculated is five. This operation is 
performed using both coefficient C{4) corresponding to both 
the last pilot symbol and coefficient C(3) corresponding to 
the last but one pilot symbol. Also in this case the 
approximation introduces a slight degradation in the overall 

20 characteristics of the system; however, the computation of 
the interpolated coefficients is on the other hand greatly 
simplified. 

The five steps required for the calculation of the five 
coefficients will be now described in a shortened way, giving 
25 the values taken by the variables, since the application 
methodology of the algorithm of Figure 3 is equivalent to the 
one described with reference to the previous cases . 

As a first step, coefficient C(7) (point C, step 1) i6 
calculated: 
30 X=4; Y=10 

Z-FLOOR [(4+10)/2]=7 

W=2-7-10«4 

Ci(7) = [Ci(4)+Cx(10) ] /2 and C Q (7) = [C Q (4) +C Q (10) ]/2 
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In the step 2, coefficient C(5) (point D,. step 2) is 
calculated: 

X^4; Y=7 

Z=FLOOR [ (4+7) /2] =5 

5 

Ci<5) = [d(3>+Cx(7) ] /2 and C Q (5) = [C Q (3 ) +C Q (7) ] /2 
In step 3, coefficient C(6) (point E, step 3) is 
calculated: 

X=5; Y=7 
10 Z=FLOOR [(5+7)/2]-6 

W=2-6-7=5 

C x (6) = [C x (5) +d (7) ] /2 and C Q (6) = [C Q (5) +C Q (7) ] /2 
In step 4, coefficient C(8) {point F, step 4) is 
calculated: 
15 X=7; Y=10 

Z=FLOOR [ (7+10) /2] =8 
W=2-8-10=6 

Ci(B) = [Ci(6)+C x (10) ] /2 and C Q (8) « [Cq(6) +C Q (10) ] /2 
In step five, coefficient C(9) (point G, step 5) is 
20 calculated: 

X-8; Y=10 

Z=FLOOR [(8+10)/2]=9 
W=2«9-10=8 

C I (9) = [C I (8)+C I (10)]/2 and C Q (9) - [C Q (8) +C Q (10) ] /2 
25 Also in this case the five coefficients C(5) C(9) 

have been obtained exclusively performing additions and 
divisions by two. 

The case in which each slot contains six pilot symbols 
(Npilot=6) and the number of channel coefficients to be 
30 calculated is equal to four, is instead illustrated in Figure 
7. Also in this case use is made of both the coefficient 
C(5) , corresponding to the last pilot symbol, and the 
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coefficient 0(4), corresponding to the last-but-one pilot 
symbol - 

The interpolation plan of Figure 7 is based on four 
consecutive steps, summarised hereinafter: 
5 As a first step, coefficient C(7) (point C, step 1) is 

calculated: 

X=5 ; Y=10 - . . . . - 

Z=FL00R [(5+10)/2]=7 
W=2-7-10=4 

10 C I (7) = [C I (4)+Ci(10)]/2 and C Q (7) - [C Q (4) +Cq (10) J /2 

In step 2, 1 coefficient C(6) (point D, step 2) is 
calculated: 

X=5; Y=7 

Z-FLOOR [ (5+7) /2] =6 
15 W=2-6-7=5 

C x (6) - [d (5) +Ci (7) ] /2 and C Q (6) - [C Q (5) +C Q (7) ] /2 

In step 3, coefficient C(8) (point E, step 3) is 
calculated: 

X=7; Y=10 
20 Z=FLOOR [(7+10)/2]~8 

W=2-8-10=6 

Ci (8) = Id (6) +d (10) ] /2 and C Q (8) « [C Q (6) +Cq (10) ] /2 
In step 4, coefficient C(9) (point F, step 4) is 
calculated: 
25 X=8; Y=10 

Z-FLOOR [(8+10) /2] =9 
W=2-9-10=8 

C X (9)-[C I (8)+C I (10) }/2 and C Q (9) - [Cq(8) +C Q (10) ] /2 
Also in this case the operations are carried out 
30 exclusively through additions and divisions by two. 

In Figure 8 there is illustrated the case in which the 
slot contains seven pilot symbols (N PIIjai =7) , said the number 
of channel coefficients to be calculated is equal to three. 
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The operation is performed in three steps using the C(6) 
coefficient corresponding to the last pilot symbol. 

As a first step, coefficient C(8) (point C, step 1) is 
calculated : 
5 X=6; Y=10 

Z=FLOOR [ (6 + 10) /2] =8 

_ _W=2-8-10 = 6 . _ . „ _ - 

Ci(8)-[C i (6)+Ci(10)]/2 and C Q (8) - [C Q (6) +C Q (10) ] /2 
In step 2, coefficient C(7) (point D, step 2) is 
10 calculated: 

X=6; Y=8 

Z=FLOOR [(6 + 8)/2]=7 
W=2-7-8 = 6 

Ci(7) = [Ci(6)+C!(8) ]/2 and C Q (7) = [C Q (6) +C Q (8) ] /2 
15 In step 3, coefficient C(9) (point E, step 3) is 

calculated: 

X=8; Y=10 

Z=FLOOR [(8+10)/2]=9 
W=2-9-10=8 

20 d (9) - [C x (8) +d (10) ] /2 and C Q <9) = [C Q (8) +C Q (10) ] /2 

In Figure 9 there is instead illustrated the last case 
in which each slot contains eight pilot symbols (N P h jOT =8) and 
the number of channel coefficients to be calculated is equal 
to two. This operation is carried out in two steps using 
25 coefficient C(6) which corresponds, in this case, to the 
last-but-one pilot symbol . 

As a first step, coefficient C(8) (point C, step 1) is 
calculated : 

X=6; Y=10 
30 Z=FLOOR [ (6 + 10) /2] =8 

W=2-8-10 = 6 

Ci(8)-[C I (6)+Ci(10)]/2 and C Q (8) = [C Q <6> +C Q (10) ) /2 
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In step 2, coefficient C(9) (point D, step 2) is 
calculated: 

X-8; Y=10 

Z^FLOOR [(8+10}/2]=9 
5 W=2-9-10 = 8 

Ci (9) - [C z (8) +Ci (10) ] /2 and C Q (9> « [C Q (8) +C Q (10) ] /2 
Also in this case the fact, of . using_ the channel 
coefficient corresponding to the last but one pilot symbol 
introduces a slight degradation in the overall system 
10 characteristics, however the computation of the interpolated 
coefficients is greatly simplified. 

The six cases previously described and illustrated with 
reference to Figures 4 to 9, can be schematised in the state 
diagram shown in Figure 10. The diagram of Figure 10, whose 
15 sequence of operations essentially depends on the value of 
the Npilot parameter, shows how it is possible to implement, 
by means of a simple state machine, a hardware or software 
device, which is the embodiment of the method of this 
invention. In the diagram there is shown only the phase 
20 component Ci(k) of the channel components, owing to the fact 
that the formulas for the computation of the corresponding 
quadrature component Cg(k) are equivalent. 

In Figure 11 there is shown instead an example of 
application, referred to the case NpmoT=6, of a second method 
25 of interpolation implemented according to the present 
invention. In the exatnple of Figure 11, for the computation 
of a plurality of channel coefficients, corresponding to the 
abscissa positions 6 to 9, use is made of the known channel 
coefficients, corresponding to the last pilot symbol, of 
30 abscissa A, of the current slot L and to the first two pilot 
symbols, of abscissas B and B+l, of the slot L+l subsequent 
to the current one. 
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To implement this second interpolation method, it is 
necessary to store an additional pilot symbol with respect to 
the cases previously described, namely the one of abscissa 11 
<B+1) , corresponding to the second pilot symbol of the 
5 subsequent slot „ 

To explain the sequence of steps illustrated in the 
application example of Figure 11, and how the same method can 
be applied also to the remaining cases, from N P i LC yr=3 to 
Npn j0T =8/ we will now analyse in detail the second 
10 interpolation method, implemented according to the present 
invention. 

This second method allows the computation, by 
interpolation, of a plurality of points comprised between a 
last channel coefficient, of abscissa A, of a current slot L, 
15 and the first two channel coefficients, of abscissas B and 
B+l, of a subsequent slot L+l . 

By way of example, use will be made within this 
algorithm of some variables that have the same meaning of the 
variables previously used with reference to the algorithm of 
20 Figure 3 : 

X: abscissa of the known left-hand point, used for the 
interpolation, a value which initially corresponds to A; 

Y: abscissa of the known right-hand point of the 
interpolation interval, a value which initially corresponds 
25 to B; 

Z: abscissa of an intermediate point between X and Y; 

W: abscissa of the right-hand extreme point, actually 
used for the interpolation corresponding from time to time 
either to point Y or to point Y+l; 
30 f (k) : ordinate corresponding to a generic abscissa k. 

The variables X and Y are initially set to be equal to 
the extremes A and B . 

Subsequently, it goes on calculating the abscissa Z, 
obviously meant as an integer number, of the intermediate 
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point between X and Y. The value of abscissa Z is rounded off 
to the higher integer through the function CEIL [(X+Y>/2]- 
Then calculation is made of the abscissa of the right-hand 
extreme point D, which is found to be equal to Y or Y+l, 
5 depending on whether the value of X + Y is even or odd, 
respectively, by means of the formula W=2«Z-X. The 
calculation is then -carried out of the ordinate f(Z) of 
intermediate point Z by arithmetic averaging between the 
ordinate of the right-hand extreme point W and the ordinate 

10 of the left-hand extreme X, which are known. 

Then the check is made on whether the point just 
calculated, of abscissa Z, corresponds to the point of 
abscissa Y-l. If this condition is met, the value of Z is 
assigned to variable X of the known left-hand point, and the 

15 operations of the preceding paragraph are applied again. In 
practice the above operations are recursively applied to each 
half- interval on the right-hand of the point previously 
calculated until the abscissa point Y-l is reached. 

Subsequently the procedure applied is by decreasing 

20 abscissas starting from the last calculated point Y-l and 
searching for a first point still to be computed. If the 
left-hand extreme A of the interval is reached, the algorithm 
is terminated, since all the points have already been 
calculated. 

25 If however a point still to be calculated has been 

found, the abscissa of which is equal to the value taken by 
variable Z, then the variable Y is assigned the value Z+l, 
which represents the first known right-hand point, and the 
search is made by subsequent attempts for the first known 

30 left-hand point of Z, a value which is assigned to variable 
X. To perform the latter operation, variable X is decreased 
each time by a unit, until a corresponding known coefficient 
is reached. 
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With the new values of X and Y, use is made again of the 
recursive algorithm for computing an intermediate point of 
abscissa Z and in cascade for the search procedure of a new 
point, still to be calculated. When all the points have been 
5 computed, the algorithm is terminated; otherwise a new 
interval X-Y is determined, to which the above recursive 
algorithm is to be applied. 

Let us now analyse how the algorithm previously 
described is applied to the interpolation plan of Figure 11. 
10 As a first step, coefficient C(8) (point C, step 1) is 

calculated: 

X=A=5; Y=B=10 

Z=CEIL [(5+10)/2]«8 

W=2-8-5=ll 

15 Ci (8) = [Ci (5) +d (11) ] /2 and Cq (8) = [C Q (5) +C Q (11) ] /2 

In step 2, coefficient C{9) (point D, step 2) is 
calculated: 

X-8; Y=10 

Z=CEIL [ (8+10) /2] =9 
20 W=2-9-8=10 

Ci(9) = CCi(8)+C x (10) ] /2 and C Q (9) = [Cq (8) +C Q (10) ] /2 

In step 3, coefficient C(7) (point E, step 3) is 
calculated: 

X=5; Y=8 
25 Z=CEIL [(5+8)/2]=7 

W=2-7-5 = 9 

Ci (7) = [Ci (5) +C t (9) ] /2 and C Q (7) « [C Q (5) +C Q (9) ] /2 
In step 4, coefficient C(6) (Point F, step 4) is 
calculated: 
30 X=5? Y=7 

Z=FLOOR [{5+7)/2]=6 
W=2-6-5 = 7 

Ci(6) = [C I (5)+C I (7) J/2 and Cq <6)«[Cq (5) +C Q (7) ]/2 
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The algorithm described in the second interpolation 
method may be regarded as a "mirror" version of the algorithm 
previously described with reference to Figure 3. As a matter 
of fact, this method, too, allows computations exclusively 
5 through additions, divisions by two and compare operations 
between registers, of any number of channel coefficients 
contained in an interval defined at the left side by a known 
coefficient (of abscissa A) and at the right side by two 
known coefficients (of abscissas B and B+l) . 

10 The two interpolation methods described above can be 

used in combination to implement a third method that might be 
advantageous in terms of speed of execution. 

In the event in which there are available as known 
extremes of the interpolation interval two channel 

15 coefficients on the left-hand, of abscissa A-l and A, and two 
known channel coefficients on the right-hand, of abscissa B 
and B+l, it is possible to apply "in parallel" both 
interpolation methods described above, in a left-hand sub- 
interval and in a right-hand sub-interval of the original 

20 interval, respectively. In fact it sufficient to calculate, 
as a first step, an average point C between the two extremes, 
A and B, of the interval, rounding up, if necessary, the 
abscissa of such point C to the next lower integer or to the 
next higher integer, indifferently, in order to define a 

25 left-hand sub-interval comprised among the abscissa points A- 
1, A and C and a right-hand sub-interval comprised among the 
abscissa points C and B, B+l. At this point it is possible to 
apply in parallel, by better exploiting therefore the 
hardware resources involved, the first interpolation method 

30 in the left-hand sub-interval and the second interpolation 
method in the right hand sub- interval, thus performing in 
parallel the computation of all the coefficients comprised in 
the interval . 
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Likewise it is possible to foresee an additional 
variation of the methods implemented according to the present 
invention; for instance it possible to apply for the 
interpolation, in place of the known channel coefficient 
immediately preceding the last known coefficient of the 
current slot or immediately following the first known 
coefficient of the _ subsequent slot, other known coefficients 
farther from the interpolation interval or a linear 
combination of the same. 

As an example of use of a linear combination of known 
coefficients the case of three pilots (N PIW>TS :=3) is analysed ♦ 
First a linear combination is performed of the three known 
pilots of the left hand group. 

CL I =C 1 (2)/2+[C I (0)+ C x (1)3/4 

CL Q =C Q (2)/2+ [Cq(0)+ 00(1)3/4, 

and the result is assigned to coefficient C(2): 
Cj(2)= Ch z 
C Q (2)= CLq 

Also Coefficient C(l) is calculated as a linear 
combination of the same pilots: 
CLx=C 3 (l)/2+[Ci(0)+ C 2 <2)]/4 
CL Q -C Q (l)/2+[C Q (0)+ Cq<2)]/4 

The result is assigned to coefficient C(l) . 
Ci(l)» CLi 
C Q (1) = CL Q 

A similar linear combination may be calculated on the 
basis of the pilots of the subsequent slot . 
CLi=C r (10)/2+[C x (ll)+ C x (12)]/4 
CL Q =C Q (10)/2+[Cg(ll)+ C Q (12)]/4 / 
and then assigned to coefficient C(10) 
Cx(10)= CI* 
C 0 (10)= CLq 
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Linear combinations are so chosen as to be implemented 
in the form of sums and shifts of bits, in a similar way as 
has been illustrated for the interpolation. 

The interpolation algorithm for the calculation of the 
5 missing coefficients can therefore be applied following the 
plan already illustrated for three pilots and reported in 
_ Figure 4 . This variation of the method is applicable- in the 
extremely frequent case of the pilots being affected by 
noise: a linear combination, which represents a weighed 
10 average of the known channel coefficients, albeit affected by 
noise, provides a more reliable estimate of the known points 
for the interpolation. 

The block diagram of Figure 12 shows a possible hardware 
implementation of a device for the channel estimation 
15 according to the present invention, in case of a transmission 
path towards a .radio base station and a spreading factor of 
the channel DPDCH equal to 256. 

The device of Figure 12, called CEU (Channel Estimation 
Unit) , allows extension of the channel estimation to channel 
20 coefficients subsequent to the pilot symbols, by implementing 
one of the interpolation methods previously described. 

The CEU unit receives at its input the channel symbols 
received after performing the operations of "descrambling" , 
»de spreading" with the code of channelling DPCCH and 
25 integration. The received pilot symbols which are supplied to 
the CEU unit are denoted by Zi(k), where subscript "i" 
indicates the signal component (i - I for the phase component 
or i = Q for the quadrature component) . 

The first operation performed by the CEU unit is the 
30 multiplication of the received pilot symbols by the reference 
pilot symbols Xq (k) . As a matter of fact, the reference 
pilot symbols are known to the receiver. Subsequently it is 
necessary to invert the sign of the phase component. The sign 
inversion of the phase component is effected through 
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multiplication by factor Mi, setting Mi = -1 for the phase 
component and Mi = +1 for the quadrature component, as shown 
in the diagram of Figure 12. 

Once the above operations have been carried out, the 
5 effect of the sequence of the pilot symbols is removed and 
the complex values obtained represent only the phase 
displacement introduced by the channel . These values are 
stored into memory 100 through the input denoted in Figure 10 
as Input__port_2 , 

10 When all the pilot symbols of the current slot and one 

or two "initial pilot symbols, of the subsequent slot, 
depending on the method to be applied, have been received and 
stored into memory 100, the process of interpolation begins 
under the control of the CLU unit 102 (CLU = Control Iiogic 

15 Unit) . 

The logic control unit 102 is a finite state machine 
(FSM) that performs a predefined sequence of operations, as a 
function of input parameter N P ilot (corresponding to the 
number of pilot symbols present in the current slot) . For 

20 instance, in case of the interpolation method previously 
described with reference to figure 3, the logic control unit 
102 carries out the sequence of operations shown in the state 
diagram of Figure 10. 

The. basic operation effected by the logic control unit 

25 102 is divided into the following three steps; 

- reading the first operand from memory 100 and loading it 
into the first register 104; 

- reading the second operand from memory 100 and loading it 
into the second register 106; 

30 - writing the arithmetic average of the two operands into 
memory 100 through its Input_jport_l • 

The computation of the arithmetic average of the two 
operands requires an addition and a division by two; the sum 
is effected in adder 108 and the division by two in block 110 
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which carries out a right-hand shift, discarding in practice 
the less significant bit (LSB) of the data resulting from 
adder 108 . 

All the calculations to be carried out only require sums 
5 and right-hand shifts, that is operations which are easily 
implement able on hardware. 

When all channel coef f icients C ± (k) , with k ranging^ from 
0 to 9, have been computed, they are sequentially read from 
the memory and supplied to the CCU unit (Channel Compensation 
10 Unit) , not shown in Figure 10, since its function is well 
known to the man skilled in the art. 

The multiplying factor shown in Figure 12 as M 2 is 
required to invert the sign of the quadrature component and 
to provide the complex conjugate of the channel estimations 
15 to the channel compensation unit (CCU) . M2 is equal to +1 for 
the phase component and to a -1 for the quadrature component. 
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